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Abstract 

By introducing a density-dependent contact term, M3Y-type interactions applicable to the 
Hartree-Fock calculations are developed. In order to view basic characters of the interactions, 
we carry out calculations on the uniform nuclear matter as well as on several doubly magic nuclei. 
It is shown that a parameter-set called M3Y-P2 describes various properties similarly well to the 
Skyrme SLy5 and/or the Gogny D1S interactions. A remarkable difference from the SLy5 and 
the D1S interactions is found in the spin-isospin properties in the nuclear matter, to which the 
one-pion-exchange potential gives a significant contribution. Affecting the single-particle energies, 
this difference may play a certain role in the new magic numbers in unstable nuclei. 

PACS numbers: 21.30.Fe, 21.60.Jz, 21.65. +f, 21.10.Dr 



I. INTRODUCTION 



Various models for nuclear structure have been developed in order to study low energy 
phenomena of the atomic nuclei. Whereas straightforward application of the bare NN 
interaction is yet limited only to light nuclei the nuclear structure seems to be well 
described by relatively simple effective interactions at low energies. Although the effective 
interactions may depend on the models, there should be basic characters in the effective 
interactions for the low energy phenomena, irrespective to the models. On the other hand, 
since the invention of the secondary beam technology, experimental data on the unstable 
nuclei have disclosed new aspects of the nuclear structure. A remarkable example is the 
dependence of magic numbers on the neutron excess 0| • In regard to the new magic numbers 
discovered near the neutron drip line, a question has been raised on a character of the effective 
interactions relating to the spin-isospin flip mode 

Mean-field theories have successfully been applied to the nuclear structure problems, 
in particular for stable nuclei. They are also useful to investigate basic characters of the 
effective interactions. However, not many effective interactions have been explored for the 
nuclear mean-field calculations so far. The Skyrme interaction Q has been popular in the 
Hartree-Fock (HF) calculations, since the zero-range form is easy to be handled. Among a 
limited number of finite-range interactions, the Gogny interaction is widely applied to the 
mean-field calculations, in which the Gaussian form is assumed for the central force. The 
parameter-sets both of the Skyrme and the Gogny interactions have been adjusted mainly 
to the data on the nuclei around the /3-stability. It is not obvious whether the available 
parameter-sets of these interactions account for the new magic numbers properly. 

In order to exploit effective interactions applicable also to unstable nuclei, guide from 
microscopic theories will be important. Brueckner's G- matrix has been a significant clue to 
studies in this course. Although microscopic approaches using the G-matrix have not yet 
been successful in reproducing the saturation properties, notable progress has been made 
recently. In the shell model approaches, microscopic effective interactions have been shown 
to reproduce observed levels remarkably well p. It should be noted, however, that the shell 
model interactions are usually specific to mass regions, and their global characters have 
not been discussed in detail, despite several exceptions (7[. The so-called Michigan 3-range 
Yukawa (M3Y) interaction Q| has been derived from the bare NN interaction, by fitting 
the Yukawa functions to the G-matrix. Represented by the sum of the Yukawa functions, 
the M3Y-type interactions will be tractable in various models. It has been shown that the 
M3Y interaction gives matrix elements similar to reliable shell model interactions [j| . More- 
over, with a certain modification, M3Y-type interactions have successfully been applied to 
nuclear reactions By using a recently developed algorithm jll|, a class of the M3Y-type 
interactions can be applied also to the mean-field calculations. Under such circumstances, 
it will be of interest to explore M3Y-type interactions and to investigate their characters 
in the mean-field framework. In this article, we shall develop M3Y-type interactions and 
investigate their characters via the HF calculations. 
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II. MODIFICATION OF M3Y INTERACTION 
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Nuclear effective Hamiltonian consists of the kinetic energy and the effective interaction, 
H = K+V- K = J2f^, V = ^. (1) 

4 Z1V1 i<j 

Here i and j are the indices of individual nucleons. It will be natural to assume the effective 
interaction Vij to be translationally invariant, except for the density dependence mentioned 
below. We consider the effective interaction having the following form, 

via = v$ + Vi2 S) + vS N) + ^i2 D) ; 

n 

(LS) = E( eSE)p TE + e SO) pTo)/ (LS) (ri2)Li2 . (si + S2); 
n 

U€ NE) Pte + 4 TNO) P T o)/f N) (n 2 ) rl 2 S 12 , 

n 

V <P») = t^) (1+x ^)p aMri)r5 (r 12 ). (2) 

The relative coordinate is denoted by ri 2 = ri — r 2 and ri 2 = [1*12 1 • Correspondingly, the 
relative momentum is defined by p 12 = (pi — p 2 )/2. L 12 is the relative orbital angular 
momentum, 

L12 = r 12 x pi 2 , (3) 
Si, s 2 are the nucleon spin operators, and Su is the tensor operator, 

S u = 4 [3(si ■ f 12 )(s 2 ■ r 12 ) - si ■ s 2 ] . (4) 

fn{f 12) represents an appropriate function of r 12 , the subscript n corresponds to the param- 
eter attached to the function {e.g. the range of the interaction), and t n is the coefficient. 
Examples of f n {r\2) are the delta, the Gauss and the Yukawa functions. P a {P T ) denotes 
the spin (isospin) exchange operator, while Pse, -Pte, -Pso and -Pro are the projection op- 
erators on the singlet-even (SE), triplet-even (TE), singlet-odd (SO) and triplet-odd (TO) 
two-particle states, respectively, which are defined by 

1 - P a 1 + P T i + Pa i-P T 

PSE-^ — , PTE-^ — , 

l-p a \-p T l + p a l + p T 

Pso-— j - . ^to-^ g - . (5) 

The nucleon density is denoted by p(r). The original M3Y interaction is represented in the 
form of Eq. (J2J), with f n {r 12 ) = e~ M " ri2 //i n r 12 and fj 2 5D ' ) = 0. As discussed in Ref. jlll ]. the 
Skyrme and the Gogny interactions are obtained by setting fn^ 12) appropriately, except for 
some parameter-sets of the Skyrme interaction in which certain terms are expressed only in 
the density-functional form. 

The saturation of density and energy is a basic property of nuclei. In developing effective 
interactions adaptable for many nuclei, it is required to reproduce the saturation property. 
However, the non-relativistic G-matrix fails to reproduce the saturation at the right density 
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and energy. Therefore, it will not be appropriate to use the G-matrix for HF calculations 
without any modification, although several HF approaches using interactions derived from 
the G-matrix were tried in earlier studies |l2j ]. The M3Y interaction was obtained so that the 
G-matrix at a certain density could be reproduced by a sum of the Yukawa functions. The 
M3Y interaction gives no saturation point within the HF theory, unless density-dependence 
is taken into account explicitly. Khoa et al. applied the M3Y interaction to nuclear reactions 
in the folding model, by making the coupling constants dependent on densities . In their 
approach the exchange terms are treated approximately. However, the exchange terms may 
contribute quite significantly to the nuclear structure. We here keep the coupling constants in 
independent of density, while introduce a density-dependent contact interaction (v in 
Eq. ©) ; as in the Skyrme and the Gogny interactions. We can then treat the exchange (i.e. 
the Fock) terms exactly with the currently available computers. It should be mentioned that 
there has been an interesting attempt to approximate the exchange terms of the interaction 
in the density-matrix expansion ^3[, although the accuracy of the density-matrix expansion 
should be checked carefully. 

We start from the Paris-potential version of the M3Y interaction This original 

parameter-set with no density- dependence is hereafter called 'M3Y-P0'. We shall modify 
this interaction so as to reproduce the saturation properties. In the isotropic uniform nuclear 
matter, matrix elements of v^ S ^ and v [™ between the HF states vanish. Therefore v^p + 
t, i2 >D ' ) determines the bulk properties such as the saturation. The range parameters for 
the Yukawa functions f^(ri 2 ) = e~^ n1 ' 12 / ' \x n r\i in are /ij" 1 = 0.25, [i^ 1 = 0.4 and 
fi^ 1 = 1.414 fm in the M3Y interaction, which correspond to the Compton wave-lengths of 
mesons with masses of about 790, 490 and 140 MeV, respectively. We do not change these 
parameters. For the longest-range part (n = 3), the coupling constants £3 , t^ E \ t^°^ and 
t^°^ are fixed to be those of the one-pion-exchange potential (OPEP), as in M3Y-P0. The 
interaction v |^ >D ' ) in Eq. (J2J) acts only on the SE and TE channels, 

W (DD) = t(DD)(1 _ x (DD) )5(ri2) p sE + t (DD) (1 + x (DD) )5(ri2) p TE (g) 

Microscopic investigations have shown that the density- dependence of the TE part is pri- 
marily responsible for the saturation ^]|, as a higher-order effect of the tensor force. While 
the interaction in the SE channel is attractive at low densities, it also has certain density- 
dependence originating in the strong short-range repulsion. Thus, a possible way of modi- 
fying the M3Y interaction may be to replace a fraction of the repulsion in the SE and TE 
channels by v^ D \ 

In addition to the saturation properties which are relevant to the central force, the LS 
splitting is significant in descri bing the shell structure of nuclei. While true origin of the 



LS splitting is not yet obvious [la], LS splittings obtained from HF calculations with the 
G-matrix interaction are too small, in comparison with the observed ones. From the HF 

(LS) 

calculations for finite nuclei, we find that v\ 2 should be about twice as strong as that of 
M3Y-P0 to reproduce the observed LS splittings. The tensor force influences the ordering 
of the single-particle (s.p.) orbits. To reproduce the observed ordering, v-^' should be 
smaller than that of M3Y-P0. We here introduce an overall enhancement factor to v\ 2 and 
an overall reduction factor to v ^ N ^ , as will be shown in Section 

In this paper we shall use two parameter-sets for modified M3Y interaction, 'M3Y-P1' 
and 'M3Y-P2', in order to show sensitivity to the parameters for some results. In M3Y-P1, 
we replace the shortest-range (n = 1) repulsive part of v±p by v^ D ^ in a simple manner. 
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We reduce both and tj™^ by a single factor, keeping the SE/TE ratio in v^ D ^ to be 
equal to t^^/tp^ in M3Y-P0, by imposing 

(TE) (SE) 

(DD) _ ? 1 1 ( 7 \ 

~ t J TE ) + t f E ) ' K ' 

The reduction factor and t( DD ^ are determined so as for the saturation density and energy in 
the nuclear matter to be typical values, as presented in the subsequent section. Characters 
of M3Y-P1 will be investigated in the nuclear matter. Although this modification is too 
simple to reproduce properties of finite nuclei, the M3Y-P1 set will be useful to clarify what 
characters arise from the original M3Y interaction, relatively insensitive to the phenomeno- 
logical modification. In the M3Y-P2 set, all t n parameters belonging to the n = 1 and 2 
channels in v l2 are shifted from those of M3Y-P0. Although we have three ranges in v\ 2 , 
the number of adjustable parameters is no greater than in the Gogny interaction, since we 

(LSI 

fix the OPEP part. We fit those parameters, together with the enhancement factor for v[ 2 
and the reduction factor for v^ N \ to the binding energies of several doubly magic nuclei. 
The resultant values of the parameters will be shown later. 



III. PROPERTIES OF NUCLEAR MATTER AT AND AROUND SATURATION 
POINT 

Basic characters of nuclear effective interactions can be discussed via properties of the 
infinite nuclear matter; in particular, properties at and around the saturation point. In 
this section we investigate characters of the M3Y-type interactions via the nuclear matter 
properties within the HF theory. In comparison, we also discuss those of the Skyrme and the 
Gogny interactions. We use the D1S parameter-set for the Gogny interaction. In most 
of the Skyrme HF approaches, the LS currents arising from the momentum-dependence of 
the central force are ignored, and the parameters are adjusted without their contribution. 
Although this treatment occasionally improves some characters of the interactions, in this 
paper we would focus on characters of the two-body interactions, rather than those of density 
functionals. For this reason we adopt the SLy5 set 18], which is devised for calculations 
including the LS currents. 

In the HF theory of the nuclear matter, the s.p. wave-functions can be taken to be the 
plane wave, 

¥w(r) = e 4k r XaXr • (8) 



Here \ a (% r ) denotes the spin (isospin) wave-function, and Q indicates the volume of the 
system, for which we will take the Q —>■ oo limit afterward. The s.p. energy for this state is 
defined as 

k 2 Q r 

e(kar) = — — + — — V / d 3 k 2 (kar, k 2 a 2 T 2 \v 12 \kar, k 2 a 2 r 2 ) . (9) 

2M (2iry Jk 2 <k FT2a2 

Energy of the nuclear matter is expressed by a function of densities depending on the spin 
and the isospin, p ra (r — p,n; a =|, [). The density variables can be converted to the total 
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density p = ^^p Ta , and the spin- and isosp in- asymmetry parameters 



err 



Ppt ~~ Ppl + P«t ~~ Pni 



P P 

Ys t p 



Vt 



J T(T 

err 



Ppt + Ppi ~~ P"T ~~ P«i 



P P 



_ _ £T _ PpT - Ppl - Pni + Pnj , . 

Vst - p p , UUj 

where a (r) in the summation takes ±1, corresponding to a =f, | (r = p,n). By assuming 
that the s.p. states are occupied up to the Fermi momentum, the density is related to the 
Fermi momentum for each spin and isospin, 

P- = ^F-- 

The total energy of nuclear matter is given by 
Q r k 2 



(27r) 3 ^ i A 1 < fcpriCT1 L 2M 

n 2 



Yj d 3 ^ d 3 k 2 (k 1 (7 1 T 1 ,li 2 a2T2\vi2\lt 1 (J 1 T 1 ,'k 2 a 2 T2) . 

(12) 



T 2 CT 2 



As already pointed out, only + ^12°^ contributes to the energy of the isotropic nuclear 
matter. In Appendix [XJ several formulae on the HF energy of the nuclear matter are derived 
for interactions expressed in the form of Eq. (j2J), with general and typical /^ c - ) (r 12 ). As well 
as for the M3Y-type interactions, the nuclear matter energies are calculated for the Skyrme 
and the Gogny interactions by using these formulae. 

In the spin-saturated symmetric nuclear matter, we have rj s = rj t = 7] st = 0, which 
indicates k^ p ] = kp P i = k-p n -\ = foxi an d PpT = Ppl = P«T = Pn[ — p/4- In this case we 
denote the Fermi momentum simply by k-p. The lowest energy for a given p normally occurs 
along this line. The saturation point is obtained by minimizing the energy per nucleon 
£ — E/A, 

f =0, (13) 

" sat. 

which yields the saturation density p (equivalently, fcp ) and energy £ . Figure [T] illustrates 
£ as a function of p for the symmetric nuclear matter with the M3Y-type as well as with 
the SLy5 and D1S effective interactions. We set M = (M p + M n )/2, where M p (M n ) is the 
measured mass of a proton (a neutron). The parameters for v^p and v^ D ^ of the M3Y-type 
interactions are listed in Table H] As mentioned above, the M3Y-P0 interaction gives no 
saturation point. We do have saturation points in M3Y-P1 and M3Y-P2 owing to v^ D \ 
Differences among the saturating forces, i.e. SLy5, D1S, M3Y-P1 and M3Y-P2, are small at 
p < p . At relatively high density (p > 0.3 fm~ 3 ), the M3Y-P1 and the M3Y-P2 interactions 
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have lower £ than SLy5 and higher than D1S. The values of /cfo and £q are tabulated in 
Table im The M3Y-P1 set has been determined so as to give A;fo — 1-36 fm and So — 16 MeV. 

In Figs. 121 and EJ contribution to £ from each of the SE, TE, SO and TO channels in 
V Y2 + t, i2 >D ' ) is shown as a function of k-p. Sum of all these channels and the kinetic energy 
(K)/A = (3/5)(fc$/2M) is equal to £ in Fig.ffl As is viewed in Fig. El the TE channel takes 
a minimum at kp = 1.3 — 1.5 fm except for M3Y-P0 and M3Y-P1, primarily responsible for 
the saturation at kpo w 1.3 fm. In the D1S interaction, the energy out of the SE channel 
monotonically goes down. This is not compatible with the presence of the strong short- 
range repulsion in the NN force, and causes an unphysical property in the neutron matter, 
as will be shown in Section IIVI Both the SO and TO channels do not contribute to £ 
significantly for p < p (> e - kp ~ k F0 ). While the SO channel becomes attractive and the 
TO channel stays small in the SLy5 and the D1S interactions, both channels are repulsive 
in the M3Y-type interactions at p > p , including M3Y-P0. A certain part of this character 
of the M3Y-type interactions comes from the OPEP part. 

The curvature at the saturation point with respect to p is proportional to the incom- 
pressibility, 



rc = k 



d 2 £ 



F dkl 



d 2 £ 



sat. 



dp 2 



sat. 



The effective mass (fc-mass) at the saturation point Mq is defined by 



<9e(kcrr) 



dk 



kpo 



(14) 



(15) 



The volume asymmetry energy corresponds to the curvature of £ with respect to r]t, 

ld 2 £ 



at 



2d V 2 



(16) 



sat. 



Analogously, the following coefficients are defined from the curvatures of £ with respect to 
i] s and r} st , 



ld 2 £ 



2drf s 



ld 2 £ 



a st = - 



sat. 



2<9r/, 



St 



(17) 



sat. 



These coefficients a s , a t and a s t are relevant to the spin and isospin responses in finite nuclei. 
In Table HTl we also compare /C, Mq, a t , a s and a st among the effective interactions. 

The incompressibility JC is sensitive to a in v^ D \ The experimental value of K, has 
been extracted from the excitation energies of the giant monopole resonances. Despite a 
certain model-dependence, most non-relativistic models are consistent with the experiments 
if K, ~ 210 MeV. For finite-range interactions, i.e. the Gogny and the M3Y-type interactions, 
a ~ 1/3 seems to give reasonable values of JC, while in the Skyrme interactions a ~ 1/6 looks 
favorable, because of the momentum-dependent terms in v l2 . The fc-mass is empirically 
known to be M * w (0.6 - 0.7)M 0. The M3Y-type interactions tend to yield slightly 
smaller Mq than the SLy5 and the D1S interactions. The volume asymmetry energy a t 
is important in reproducing global trend of the binding energies for the Z ^ N nuclei. 
From empirical viewpoints a t ~ 30 MeV seems appropriate, as is fulfilled in the M3Y-type 
interactions under consideration. 

The a s and a s t coefficients are relevant to the spin degrees of freedom. The kinetic energy 
has a certain contribution to a s and a st , as well as to a t , which amounts to about 12 MeV at 
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p ~ po equally for a t , a s and a st . The interaction v& + f 12 D ^ gives rise to the rest of these 
coefficients. Both the M3Y-type interactions have similar tendency with respect to these 
coefficients. It is remarkable that a st is substantially larger in the M3Y-type interactions 
than a s . As is suggested by close a s and a st values between M3Y-P1 and M3Y-P2, the 
original M3Y interaction already carries this feature. In particular, the OPEP part included 
in the M3Y-type interactions plays a significant role, increasing a st by about HMeV. On 
the other hand, a s and a s t are comparable in the Gogny D1S interaction, and we have even 
a s > a st in the Skyrme SLy5 interaction. In the SLy5 case, a st is close to the value due only 
to the kinetic energy. 

Global characters of the spin and isospin responses are customarily discussed in terms of 
the Landau parameters. Formulae on the Landau parameters at the zero temperature are 
given in Appendix El We compute the parameters of Eq. (jB22|) . The results are shown in 
Table ITTT1 It is remarked that the M3Y-P1 and M3Y-P2 interactions give similar results. The 
Qt and the g' e parameters are closely related to the a s and the a st coefficients, respectively. 
It has been known that go is small, while g' Q should be relatively large |20j . Although it is 
not easy to extract precise values of the Landau parameters from experimental data because 
they could depend on the interaction forms, qualitative trend will not depend on effective 
interactions. The M3Y-type interactions seem to have reasonable characters on the spin 
and isospin responses, while SLy5 and D1S do not, although the spin and isospin natures 



of the Skyrme interactions seem to be improved if the LS currents are ignored |21|. It is 
likely that the difference in these coefficients may significantly influence predictions of the 
spin and isospin responses of finite nuclei. 

IV. PROPERTIES OF ASYMMETRIC NUCLEAR MATTER AND NEUTRON 
MATTER 

We turn to the asymmetric nuclear matter. In Fig. HJ energies per nucleon £ are depicted 
as a functions of p for the spin-saturated (i.e. r) s = r] st = 0) nuclear matter with rj t = —0.2 
and —0.5. The results from the M3Y-type interactions are compared with those of the 
Skyrme and the Gogny interactions. Energies of the spin-saturated neutron matter (i.e. 
rjt = —1) are presented in Fig. Results from a microscopic calculation in Ref. [22| are 
also shown as a reference. Although the dependence on the interactions is not strong at low 
densities even for the neutron matter, it becomes stronger at p > 0.2 fm as \r] t | increases. In 
the D1S result for the neutron matter, £ has a maximum at p ~ 0.6 fm and goes to —00 as 
p — > 00. This unphysical behavior arises from 

x (dd) = ! in 

the D1S set, which implies no 

density-dependence in the SE channel (see Eq. (JBJ)). This could also give rise to a problem 
in practical calculations for finite nuclei. With the SLy5 interaction £ goes up rapidly at 
any r] t , because of the momentum- dependence of the interaction. In contrast to them, the 
M3Y-type interactions give moderate £ for the neutron matter. The microscopic energy of 
Ref. |22J lies between those of M3Y-P1 and M3Y-P2. It will be possible, if necessary, to 
adjust the parameters of the M3Y-type interactions to the microscopic results. 

V. PROPERTIES OF DOUBLY MAGIC NUCLEI 

We next discuss properties of doubly magic nuclei in the HF app roximation. In the 



calculations for finite nuclei, we use the algorithm presented in Ref. where the following 
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s.p. bases are employed, 



Raljir) = M al3 r e+2 ^ exp[-(r/z/ a ) 2 ] . (18) 

expresses the spherical harmonics. We drop the isospin index without confusion. 
The index a indicates p a (a non-negative integer) and u a , simultaneously. By choosing p a 
and v a appropriately, these bases span the space equivalent to that of the harmonic-oscillator 
(HO) bases, as well as they can form the Kamimura-Gauss (KG) basis-set j^]. Without 
parameters specific to mass number or nuclide such as hu, a single set of the KG bases 
is applicable to wide range of nuclides. In the following calculations we apply the hybrid 
basis-set |ll( for the nuclei with A < 50, in which an HO basis is added to the KG basis-set, 
while the HO basis-set with N osc < 15 and Tuuj = 41.2A~ 1 / 3 MeV for heavier nuclei. 

In finite nuclei the non-central forces are important as well. In the M3Y interaction, 
the LS force v^f 1 ^ and the tensor force fi™"* are taken by setting /i, LS ' ) (r 12 ) = e~ Mnri2 //i n r 12 
and /^ TN ^(r i2 ) = e _/inri2 //i n r 12 in Eq. (j2J). We here fix the range parameters as in v^p; 
/if 1 = 0.25 fm, jj,^ 1 = 0.4fm for v^ S \ and /if 1 = 0.4fm, jj^ 1 = 0.7 fin for v^ N \ The 
coupling constants in the M3Y-P2 set are tabulated in Table I1V1 together with those in the 

(IjS) 

original M3Y-P0 set. In M3Y-P2, the enhancement factor for v l2 is taken to be 1.8 and the 
reduction factor for v\ 2 to be 0.12. The binding energies and the rms matter radii obtained 
from the HF calculations with M3Y-P2 are shown in Table 13 in comparison with those of 
the SLy5 and the D1S interactions, as well as with the experimental data. The one-body 
terms of the center-of-mass (cm.) energy are removed before iteration. The contribution of 
the two-body terms is subtracted from the convergent HF wave-functions, in the D1S and 
the M3Y-P2 results. There are also spurious cm. effects in the matter radii, 

(r 2 ) = ^E((^-R) 2 ) = jU^)-(R 2 ) 

(19) 

The first term in the right-hand side is expressed by one-body operators with a correction 
factor (1 — 1/A). We need two-body operators for the second term. For the D1S and the 
M3Y-P2 interactions we fully remove the cm. contribution according to Eq. ()19j) . For the 
SLy5 interaction we use only the one-body terms with the correction factor, ignoring the 
two-body terms in Eq. (|T9"]). as in calculating the energies. 

Wave-functions of the doubly magic nuclei are considered to be well approximated in 
the spherical HF approaches. It should still be noted that correlations due to the residual 
interaction could influence their properties. Therefore we do not pursue fine tuning of the 
parameters. As shown in Table the M3Y-P2 set is fixed so as to reproduce the measured 
binding energies of the doubly magic nuclei, including 90 Zr, within about 5MeV accuracy. 
The binding energies of these nuclei obtained from the SLy5 and the D1S interactions are 
in agreement with the experimental data within 3 MeV, slightly better than M3Y-P2. We 
do not have to take this difference seriously, before evaluating influence of the residual 
interactions. As well as the binding energies, the rms matter radii of these nuclei are 
reproduced by the M3Y-P2 set similarly well to the other available interactions. In Table IVT1 
we present the neutron s.p. energies e n (0p 3 / 2 ) and e n (0pi/ 2 ) around 16 0. The enhancement 
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factor for v[ 2 in the M3Y-P2 set has been adjusted approximately to the experimental 
value of this s.p. energy difference. The reduction factor for t>j™' ) has been determined so 
as to reproduce the s.p. energy ordering for 208 Pb. Without this reduction factor, the orbits 
with higher I have too high energies. The resultant s.p. levels in 208 Pb with M3Y-P2 are 
depicted in Fig. |H1 The levels obtained from D1S and the experimental s.p. levels are also 
shown. The overall level spacings are related to Mq shown in Table |HI In the usual HF 
calculations the level spacings tend to be larger than the observed ones, and it is not (and 
should not be) remedied until the correlations due to the residual interaction (or the cu-mass) 
are taken into account (19:]. This is also true in the present case. We find that M3Y-P2 
yields as plausible s.p. levels as D1S does. We thus confirm that the M3Y-P2 interaction 
well describes global nature of stable nuclei. 



VI. SINGLE PARTICLE LEVELS IN N = 16 ISOTONES 

In the preceding section we have shown that the M3Y-P2 interaction reproduces the prop- 
erties of the doubly magic nuclei to a similar accuracy to the SLy5 and the D1S interactions. 
At a glance, the spin-isospin characters in the nuclear matter, which have been discussed 
in Section IIHI via a st and g' e , do not seem to influence the nuclear properties around the 
ground states. However, the spin and isospin characters influence s.p. energies of finite nu- 
clei. Thereby they may affect even the ground state properties. In this section we illustrate 
this point by the neutron orbits in the N = 16 isotones, following the arguments in Ref. 0], 
although precise studies in this line are beyond scope of this paper. 

As was suggested in Ref. the proton-number (Z) dependence of the neutron s.p. 
energy e n {Qdz/2) relative to e n (lsi/2) can sizably be affected by effective interactions. Figured 
depicts Ae n = e n (0d3/2)— e n 0- s i/2) obtained from the spherical HF calculations in the N = 16 
isotones. Though it is not obvious whether the ground states of all of these isotones are well 
approximated by the spherical HF wave-functions, it is meaningful to see the s.p. energies, 
which often give an indication to magic or submagic numbers. For D1S we reduce the number 
of bases in Eq. (|18|) to avoid instability occurring for some N = 16 nuclei, which probably 
relates to the unphysical behavior in the neutron matter. It is found that, if viewed as a 
function of Z, Ae n strikingly depends on the interactions. With the M3Y-P2 interaction, 
Ae n increases as Z goes from Z = 14 to Z = 8. We have confirmed |28[ that even M3Y-P1 
(with appropriate v^f^ and ) shows similar behavior and that a significant part of this 

(V) 

feature originates in the OPEP part in v\ 2 . It is thus suggested that this behavior of Ae n 
is correlated to the spin-isospin property in the nuclear matter. 

For comparison, we also show the s.p. energies obtained from the reliable shell model 
interaction for the sc?-shell nuclei, the so-called USD interaction For this purpose we 
define effective values of the s.p. energies for each nucleus A from the shell model space and 
interaction, which correspond to those of the spherical HF calculations, as 

e v n SD (r, a) = 6^0; 17 o) + E(iv» A (2j+ 2 1 J ) ^ + 1) to"; ^ USD b"; J) , m 

where the sum with respect to j' runs over the valence orbits. For (Nj)a, we assume that 
the nucleons occupy the s.p. orbits from the bottom, according to e(j). From these s.p. 
energies we obtain ^ie^ SD = e^ SD (0d 3 /2', A) — e^ SD (lsi/ 2 ; A) for individual nucleus. This 
definition is equivalent to the effective s.p. energies in Ref. |3| for the Z < N(— 16) nuclei. 
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The Z\e^ SD values are also shown in Fig. [7| It is noted that in the shell model approaches 
the nucleus-dependence of the s.p. wave-functions is not fully taken into account. Effects 
of rearrangement in the wave-functions of the deeply bound orbits are renormalized into 
the interactions among the valence nucleons. In contrast, in the HF approaches the s.p. 
wave-functions are determined self-consistently, from nucleus to nucleus. Therefore the shell 
model s.p. energies do not agree with their HF counterparts. However, there should be 
qualitative correspondence, which arises from basic characters of the effective interactions. 
It is remarked that the M3Y-P2 interaction has the same trend of Ae n , in terms of the 
Z-dependence, as the USD interaction. It has been suggested |3] that the interaction in 
the (<r • er) (r ■ r) channel, which will be linked to a s t or to g' E , is significant to the magic 
numbers in highly neutron-rich nuclei, and that the Z-dependence of the s.p. energies in 
this region could be relevant to the new magic number N = 16 [2j|. The present results are 
fully consistent with the arguments in Ref. |3[, although we cannot draw conclusions on the 
magic number problem without assessing influence of the residual interactions. 



VII. SUMMARY AND OUTLOOK 

We have developed effective interactions to describe low energy phenomena of nuclei. 
Starting from the M3Y interaction, we introduce a density-dependent contact term and 
modify several parameters in a phenomenological manner, whereas maintaining the OPEP 
part in the central force. In order to view basic characters of the interactions, the Hartree- 
Fock calculations are implemented for the infinite nuclear matter (for which useful formulae 
are newly derived) and for several doubly magic nuclei. We have shown that a parameter-set 
called M3Y-P2 describes their properties plausibly. The properties which are well treated 
by the Skyrme SLy5 and/or the Gogny D1S interactions are also reproduced by the M3Y- 
P2 interaction. However, a remarkable difference is found in the properties relevant to the 
spin degrees of freedom in the nuclear matter. The M3Y-type interactions seem to give 

(C) 

reasonable spin and isospin properties, in which the OPEP part contained in v\ 2 plays a 
significant role. We have also shown that the difference in the spin-isospin property affects 
the s.p. energies in finite nuclei to a considerable extent. It will be interesting to apply 
extensively the M3Y-type interactions, particularly to the magic number problems far from 
the /3-stability. 

Although the M3Y-P2 interaction seems to have various desired characters, there still 
remains a certain room for further tuning of the parameters. It should be noted that this 
parameter-set will not be a unique choice to reproduce the properties of the nuclear matter 
and the doubly magic nuclei. Effective interaction might not be constrained sufficiently 
only from the HF calculations. The pairing effects in nuclei give valuable information on 
the effective interaction, primarily on the SE channel. Comparison of the matrix elements 
with reliable shell model interactions will also be helpful, if the core polarization effects are 
treated appropriately. These points will be discussed in future publications. 
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TABLE I: Parameters of central forces (including v l2 ) in the original and modified M3Y inter- 
actions. See text for the \x n parameters. 



parameters 




M3Y-P0 


M3Y-P1 


M3Y-P2 


t m 


(MeV) 


11466. 


8599.5 


8027. 


t (TE) 


(MeV) 


13967. 


10475.25 


6080. 


t |so) 


(MeV) 


-1418. 


-1418. 


-11900. 


4 TO) 


(MeV) 


11345. 


11345. 


3800. 


4 SE) 


(MeV) 


-3556. 


-3556. 


-2880. 


4 TE) 

2 


(MeV) 


-4594. 


-4594. 


-4266. 




(MeV) 


950. 


950. 


2730. 


.(TO) 


(MeV) 


-1900. 


-1900. 


-780. 


(SE) 
r 3 


(MeV) 


-10.463 


-10.463 


-10.463 


(TE) 
r 3 


(MeV) 


-10.463 


-10.463 


-10.463 


(SO) 
E 3 


(MeV) 


31.389 


31.389 


31.389 


(TO) 
l 3 


(MeV) 


3.488 


3.488 


3.488 


a 






1/3 


1/3 


t (DD) 


(MeVfm) 


0. 


1212. 


1320. 


X (DD) 






0.09834 


0.72576 



TABLE II: Nuclear matter properties at the saturation point. 







M3Y-P1 


M3Y-P2 


SLy5 


D1S 




(fm) 


1.358 


1.340 


1.334 


1.342 


£o 


(MeV) 


-15.99 


-16.14 


-15.98 


-16.01 


K 


(MeV) 


225.7 


220.4 


229.9 


202.9 


M$/M 




0.641 


0.652 


0.697 


0.697 


a t 


(MeV) 


30.35 


30.61 


32.03 


31.12 


a s 


(MeV) 


20.81 


21.19 


37.47 


26.18 


a st 


(MeV) 


37.63 


38.19 


15.15 


29.13 
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TABLE III: Landau parameters at the saturation point. 





M3Y-P1 


M3Y-P2 


SLy5 


D1S 


fo 


-0.370 


-0.357 


-0.276 


-0.369 


fl 


-1.078 


-1.044 


-0.909 


-0.909 


h 


-0.381 


-0.436 


0.0 


-0.558 


h 


-0.191 


-0.210 


0.0 


-0.157 


fo 


0.525 


0.607 


0.815 


0.743 




0.537 


0.635 


-0.387 


0.470 


A 


0.250 


0.245 


0.0 


0.342 


fs 


0.101 


0.096 


0.0 


0.100 


9o 


0.046 


0.113 


1.123 


0.466 


9i 


0.372 


0.273 


0.253 


-0.184 


92 


0.199 


0.162 


0.0 


0.245 


93 


0.088 


0.078 


0.0 


0.091 


9'o 


0.891 


1.006 


-0.141 


0.631 


9'i 


0.230 


0.202 


1.043 


0.610 


9'2 


0.073 


0.040 


0.0 


-0.038 


9's 


0.008 


-0.002 


0.0 


-0.036 



TABLE IV: Parameters of non-central forces in the original and modified M3Y interactions, 
text for the fj, n parameters. 



parameters M3Y-P0 M3Y-P2 

4 LSE) (MeV) -5101. -9181.8 

j(LSO) ( MeV ^ _ 1897 _ -3414.6 

4 LSE) (MeV) -337. -606.6 

4 LSO) (MeV) -632. -1137.6 

tf NE) (MeV) -1096. -131.52 

4 TNO) (MeV) 244. 29.28 

4™ E) (MeV) -30.9 -3.708 

4™ 0) (MeV) 15.6 1.872 
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TABLE V: Binding energies and rms matter radii of several doubly magic nuclei. Experimental 
data are taken from Refs. 24, 25|, 2f|. 









Exp. 


M3Y-P2 


SLy5 


D1S 


16Q 


—E 


(MeV) 


127.6 


127.1 


128.6 


129.5 




vV> 


(fm) 


2.61 


2.60 


2.59 


2.59 


40 Ca 


-E 


(MeV) 


342.1 


338.7 


344.3 


344.5 




V(r' 2 ) 


(fm) 


3.47 


3.37 


3.29 


3.36 


48 Ca 


-E 


(MeV) 


416.0 


411.8 


416.0 


416.8 




y/(T*) 


(fm) 


3.57 


3.52 


3.44 


3.50 


90 Zr 


-E 


(MeV) 


783.9 


778.7 


782.4 


784.5 






(fm) 


4.32 


4.25 


4.22 


4.23 


132 Sn 


-E 


(MeV) 


1102.9 


1098.1 


1103.5 


1102.9 






(fm) 




4.79 


4.77 


4.76 


208p b 


-E 


(MeV) 


1636.4 


1635.8 


1635.2 


1638.1 




V(r 2 ) 


(fm) 


5.49 


5.53 


5.52 


5.51 



TABLE VI: LS splitting around 16 0. Experimental data are extracted from Refs. 24, 





Exp. M3Y-P2 


SLy5 


D1S 


e n (0p 3/2 ) (MeV) 


-21.8 -22.6 


-20.6 


-22.3 


e n (0 Pl/2 ) (MeV) 


-15.7 -16.2 


-14.4 


-15.9 
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FIG. 1: Energies per nucleon £ = E/A in the symmetric nuclear matter for several effective 
interactions. The thick dotted, dot-dashed and solid lines represent the results with the M3Y-P0, 
M3Y-P1 and M3Y-P2 interactions, respectively, while the thin dashed and solid lines those with 
the SLy5 and D1S interactions. 
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150 -i 




"I 1 1 1 1 1 
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k F [fm" 1 ] 



FIG. 2: Contribution of the SE and TE channels to 8. See Fig. ^ for conventions. 
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TO channel 



i 1 1 1 1 — 
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FIG. 3: Contribution of the SO and TO channels to 8. In both channels, the results of M3Y-P0 
is equal to those of M3Y-P1, which are presented by the dot-dashed line. See Fig. ^ for the other 
conventions. 
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FIG. 4: Energies per nucleoli £ = E/A in the asymmetric nuclear matters with r/t = —0.2 and 
—0.5 for several effective interactions. See Fig. [I\for conventions. 
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FIG. 6: Single-particle energies for 208 Pb. Experimental values are extracted from Refs. |24l. I27 1 
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Z 

FIG. 7: Ae n for the N = 16 isotones. The thick solid, dotted, thin solid and dashed lines correspond 
to the results with the M3Y-P2, USD, D1S and SLy5 interactions, respectively. 
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APPENDIX A: ANALYTIC FORMULAE FOR NUCLEAR MATTER ENERGY 



In this Appendix we derive formulae concerning the interaction part of Eq. (jl2j) . The 
form of Eq. © is assumed for i> 12 . 

Each term of v^p is expressed as f^\ri 2 )OcrO T . Its non-antisymmetrized matrix element 
in the plane wave states of Eq. (jSJ) is evaluated as 

= ffj dVAe^-KK^-H)^ /(0(r 12 ) {a' x a' 2 \O a \a x a 2 ) (t[t' 2 \0 » T |nr a > 

= rf 3 i?rf 3 r 12 e^ K - K ') R+ ^ k -- k '-)- r -/f)(r 12 ) (o^^ki**) (r^iO.inr.) 

= ^ <*k,k< /f } (|ki 2 - k' 12 |) Ka 2 |aki^ 2 ) (r^|a|rir 2 ) , (Al) 

where R = (r x + r 2 )/2, r 12 = r x - r 2 , K = k x + k 2 , K' = k' x + k 2 , k i2 = (kj - k 2 )/2, 
k' 12 = (k' x — k' 2 )/2 and /(g) is the Fourier transform of f(r), 

/(?) = Jd 3 rf(r)e- l *\ (A2) 

The density- dependent interaction t> is also handled in a similar manner, since the density 
behaves like a constant in the nuclear matter. For the Hartree term we have (kicriTi) = 
(k^cr^) and (k 2 a 2 r 2 ) = (k 2 cr 2 r 2 ), while (kiCTiTi) = (k 2 cx 2 T 2 ) and (k 2 cr 2 r 2 ) = (k^cr^) for the 
Fock term. Therefore both terms satisfy K = K'. For the relative momentum the Hartree 
term (the Fock term) yields k 12 — k' 12 = (k 12 — k' 12 = 2k 12 ). Contribution of the two-body 
interaction to the nuclear matter energy is obtained by integrating / in Eq. (jAl|) up to the 
Fermi momenta. 

We here consider general cases where the Fermi momentum may depend on spin and 
isospin. In order to take into account the spin-isospin dependence, we integrate / in the 
range k x < k-p\ and k 2 < k F2 . The integration is immediately carried out for the Hartree 
term, as far as f(r i 2 ) is momentum- independent, since the integrand depends neither on ki 
nor on k 2 , 

W H (fc F1 , k F2 ) = / d 3 h / d 3 k 2 f(0) = -—k 3 F1 k 3 2 f(0) . (A3) 

Jki<kpi Jk2<kp2 y 

For the Fock term contribution, the integral with respect to ki and k 2 is converted to the 
one with respect to K and ki 2 . We here assume /cfi < kp 2 without loss of generality, owing 
to the symmetry W(kpi, k^) = W(A:f2,*fi)- Handling the range of integral carefully, we 
obtain the following expression, 



W F (A; F1 ,A; F2 ) = / d 3 h f d 3 k 2 f(2k 12 ) 

jfcl<fcpi Jfc2<fcp2 

/ dk 12 -k 3 F1 k 2 12 f(2k 12 ) 

JO 6 

/ dk 12 --(k 2 F2 - k 2 F1 ) 2 k 12 + -(k F1 + k 3 2 )k 

J(k F2 -k F1 )/2 12 6 

-i(k F1 + k F2 )k 3 12 + ~k 5 12 } 7(2*12 



(A4) 
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These formulae are general to multi-component uniform Fermi liquids with equal masses. 
In handling the spin-isospin degrees of freedom, we rewrite the central force in Eq. (J2J) as 



u 12 



E(4 W) + t^P. - tWp T - t^P a P T )fP(r 12 ) . (A5) 



The relations between the coupling constants are 

,(SE) _ ,(W) _ ,(B) _ .(H) , f (M) ,(TE) _ ,(W) , ,(B) , .(H) , ,(M) 
*(SO) _ f (W) _ f (B) , ,(H) _ ,(M) .(TO) _ .(W) , .(B) _ ,(H) _ .(M) / A fi N 

After summing over the spin-isospin degrees of freedom, the interaction energy is given by 

\ / n/r) \q Z_^i L^i [\ n ' b n KJ cr 1 (j 2 h n u t\t 2 h n u o- 1 <T2 u t 1 t 2 J ' v n V^VTyai i <^¥t 2 <j 2 J 

Z l Z7r J n a 1 cr 2 T 1 T 2 

, f f (M) , f (tt) x _/( B )A _/CW)j; x )\V F (kv K, ^1 (A7) 

' \ L n ' L n v cr\u 2 L n u tit 2 L n u u 1 a 2 ^>T 1 T 2 J ' *n \ h, F~nai J "'FT2>T2j\ ■\- n -'J 

In Eq. flA7|) we regard the sum over n to include v^ D \ It is noted that A = Qp, which is 
used to obtain the energy per nucleon S. 

We next calculate the W functions for typical interaction forms. 

1. 5 interaction 

If fij '12) = <K r i2), f(l) = 1 an d therefore we have 

W H (A; 1 , A*) = W F (fci, fca) = ^-kfk 3 , . (A8) 

2. p- dependent 6 interaction 

Since the density is a constant in the uniform nuclear matter, the W functions for 
fij '12) = p a ^(i"i2) are similar to the above case, 

W H (A; 1 , fc 2 ) = W F (h, k 2 ) = -^p a k\kl . (A9) 

Note that p is a function of the Fermi momenta, when we take derivatives of the W 
functions. 

3. Gauss interaction 

For /(ri 2 ) = e"^ ri2 ) , we have /(g) = (y^rr/jj,) 3 e'^ 2 ^ , deriving 

W^k u k 2 ) = 1 -^(^-)\fkl, (A10) 



and 

W F (h,k 2 ) = ^# 



■7 r 



p |(& 2 - k x k 2 + k 2 2 - 2p 2 ) e ~(^) 

-(hi + A;iA; 2 + A; 2 - 2p 2 ) e'^^^ 
-(*? + kl) Erfc(^i) + (kl - k\) Erfc(^i) + y/jk* 



(All) 
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where 

Erfc(ar) = / e~ z 'dz. (A12) 



In Eq. flAll|) we have postulated fci < k 2 again. 



4. Yukawa interaction 

For the Yukawa interaction we set /(Vi 2 ) = e~^ Tl2 / /iri 2 , leading to /(g) = 4tc/ fi(fi 2 + 
q 2 ). This yields 

64tt 3 

"V 

and 



W H (fc 1 ,A; 2 ) = ^^ 2 3 , (A13) 



W F (h,k 2 ) = ^[Ak 1 k 2 {3(k 2 1 + kl)-ii 2 } 



-16// | (A; 3 + arctan ^ 1 * — j - (A| - A; 3 ) arctan ^ 2 - kl 
- {3(k 2 2 - kl) 2 - 6fi 2 (kl + k 2 2 ) - /i 4 } In 



I + (** + ^) 2 " fAl4 N 

fi 2 + (k 2 -k 1 y ' u 



5. Momentum- dependent 6 interaction 

In the Skyrme interaction we have momentum-dependent terms with the form 
|{Pi2^( ri2 ) + < K r i2)Pi2} an d P12 • ^( r i2)Pi2- The former operates only on the even 
channels and yields 

4vr 2 

W H (fc 1; k 2 ) = W F (h, k 2 ) = —k 3 k 3 (k 2 + k\) . (A15) 

15 

The latter acts on the odd channels, giving 

4tt 2 

W u (h, k 2 ) = -W F (fc 1; k 2 ) = —k\kl{kl + k\) . (A16) 

15 

The incompressibility K and the spin-isospin curvatures a t , a s , a st are expressed by the 
derivatives of the W functions. 

The single-particle energy e(kcrr) defined in Eq. (J9J) is also expressed by the derivative of 
the W functions. We first rewrite the integral in Eq. (|12|) as 

/ d 3 k[ / d 3 k 2 {k' 1 a 1 T 1 ,k 2 a 2 T 2 \v 12 \k' 1 aiTi,k 2 a 2 r 2 ) 

= An [ k' 2 dk[ [ d 3 k 2 {k' 1 a 1 T 1 ,'k 2 a 2 T 2 \vi 2 \k' 1 a 1 T 1 ,k 2 a 2 T 2 } . (A17) 

JO Jko<k w ^„„„ 



This immediately gives 




/ d 3 k[ / d 3 k 2 (k' 1 aiT 1 , k 2 o 2 T 2 \vi 2 \k' 1 a 1 n,k 2 a 2 r 2 ) 

OK\ Jk[<ki J k2<kYT 2 a2 

= Airkl d 3 k 2 (k' 1 o-iTi,k 2 a 2 T 2 \vi 2 \k' 1 aiTi,k 2 a 2 T 2 ) . (A18) 

Jk 2 <k FT2tT2 



24 



Therefore, 
e(k 1 cr 1 n) 



k 2 1 1 



2M (27r) 3 4vrA; 2 ^ 

Y+(W) , f (B)r _+( h )a -f( M )x A ^W^fri 

+ (if) + tW<W 2 - t^6 TlT2 - t^5 aia2 5 T1T2 ) d^ih, k F 



where we use the short-hand notation 



_d_ 



(A19) 



(A20) 



It is now obvious that the effective mass of Eq. (fTHjl is expressed by using the second 
derivative of the W functions. 



APPENDIX B: LANDAU PARAMETERS FOR SYMMETRIC NUCLEAR MAT- 
TER 

Let us denote the occupation probability of the s.p. states of Eq. (jSJ) by n T(T (k). The 
nuclear matter energy of Eq. (jA7|) can be rewritten as 



(V) (V)n + (V) } 



(V)h 1 



(Bl) 



n 2(2tt) 6 ^ ^ E n n<xi( k i)^r 2CT2 (k 2 )/n(0) 



n u\(72T\T2 kik2 

f t (w) +t (B) 5 _m§ -tW>6 5 ) (B2) 
= 9(9v\e n ri<n (ki)r2 T2(7 2 (k 2 ) fn{2k 12 ) 

\l Z{Z7T) n ai(72TlT2 klk2 

. ( f (M) , f (H)r _f( B )A -f^S 8 \ (Bl) 

The Landau coefficient is defined by 

F^ a2 (h, k 2 ) = 2 ^f 1 d(k ■ k 2 ) P e (k ■ k 2 ) - . (B4) 

1 1,2 2 2 J-i 5ra TlCTl (ki)5ra T2(72 (k 2 ) 

For the interaction independent of momentum and of density, it is straightforward to write 
down the coefficients of Eq. ()B4j) in terms of /, within the HF theory at the zero temperature. 
Noticing that p also depends on n Tcr (k), we evaluate contribution of the density-dependent 
5 interaction (1 + P a )p a S(r 12 ) to F^ ai>T2a2 {k u k 2 ) as 

[^pV 2 {p 2 - E & + * (DD) (e pI - E pi) } 

+ap a ~ 1 [2p - p Tiai - p T2f72 + x (DD) (p ai + p a2 - p Tl - p T2 )} 

+p a {l - 5 T1T2 5 aia2 + x^ D \5 ai(72 - 5 T1T2 )}} , (B5) 
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where p a = J2 T Pra and p T = J2a Pra- Apart from the spin and isospin degrees of freedom, the 
momentum-dependent 5 interactions |{Pi2<K r i2) + ( K r i2)Pi2} an d P12 • <K r i2)Pi2 contribute 
to F^^kuh) by 

1 ^S±M-^V (B6) 



(2vr; 

In characterizing effective interactions, we view the Landau coefficients for the symmetric 
nuclear matter, where p rcr = p/4 for any r and a. While formulae for the Landau parameters 
were derived for the Skyrme interaction in Ref . 21 1 and for the Gogny interaction in Ref . |3(| , 
we here derive expressions for interactions with the form of Eq. (J2J) in more general manner. 
It is customary to transform the (r, a) variables into the following ones, 

1 ■ • ■ p T +P i +n T +n I , 

t ■ ■ ■ P t +P I -n t -n i , ^ 
s - ■ • p | — p i +n t — n i , 
st ■ ■ ■ p t —p I — n t +n i . 

Since ° — H T T = Z)cr( crr ) = J2 T (< JT ) = 0, & U the off-diagonal coefficients with respect to 
(l,t, s, st) vanish. The diagonal coefficients are redefined as 





1 

" 16 


E 

(J1(T2T1T2 




F t {£ \h,k 2 ) 


1 

" 16 


E 




F^XhM) 


1 

~ 16 


E 


a^F^^kz), 


F { J\k u k 2 ) 


1 

" 16 


E 

<T1<T2T1T2 


a 1 T 1 a 2 T 2 F^l UT2a2 (kt, k 2 



(B8) 

The Hartree terms of the momentum- and density-independent interactions yield 

E(4t (W) +2t (B)_ 2t (H)_ e))/n(0); 
n 

E( _ 2t (H)_ t (M) )/n(0); 

n 

, 6 E(24 B) -e ) )/n(o), 

n 

E(-4 M) )/n(0), (B9) 

while the Fock terms 



F[%{hM) 


Seo 
4(2tt) 


F t %(h,k 2 ) 


4(2vr) 


F^%(k h k 2 ) 


4(2tt) 




Seo 
4(2tt) 



6 
n 



F^k 2 ) = ^E(< M) + 24 H) -2e)-e))Gf(^,^), 



4(2vr) b n 
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F^Ah, k 2 ) = E(-4 W) ) k 2 ) , (BIO) 

where 

Of 4. 1 rl 

G®(h, k 2 ) = — ^— jf_ i d(ki • k 2 ) P,(ki • k 2 ) f n (2k 12 ) . (Bll) 
Contribution of the density-dependent interaction t^ DD ^(l + x^ DD ^ P a )p a 6 (i"i 2 ) is given by 

,3(a + l)(a + 2) ^ Q 





<^0 , 

" 4(2tt) 6 


F t %»(h,k 2 ) 


~ 4(2tt) 6 
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~ 4(2tt) 6 




5m 
4(2tt) ( 



■t (DD) p a . (B12) 

For momentum-independent interactions such as the Gogny interaction and the M3Y- 
type interactions, the Landau coefficients are obtained by Fi\ki,k 2 ) = F^(ki,k 2 ) + 
Fyl(ki, k 2 ) + ^i dd(^1j ^2), and so forth. The momentum-dependent 5 interactions yield 

pW „ , . 1 ( x kl + kl . . . \ |3tS MD) 

-^l,MDV"'l) k 2 ) - ^ 2n y \ £0 2 ^l fc l fc 2 I X I 5t (MD) > 

H" , , 1 ^ fc i 2 + A; 2 , , , \ fti MD) (-2xS MD) -l) 

W*i, **) - 2— ~ 6ilhh ) X \ 4 MD )(2x 2 MD) + 1) ' 



kl 


+ kl 




2 


kl 


+ kl 




2 


kl 


+ A; 2 2 




2 


kl 


+ A; 2 2 


2 



zpW „ ^ 1 (, kl + ki . , . \ /4 MD) (2xi MD) -l) 

- 3(2^6^*— a" -fe*i*>J X \4 MD )(2x 2 MD) + l) ' 

**) = ^ (^^±M - tofc,) x I (^ m) ) , (B13) 

where the upper row corresponds to the even channel interaction |ij ^(l + 
^(MD)p^^|p 2 ^^^ r ^^-j _|_ ^^ ri2 ^p2 2 j^ w ] 1 ji e l OW er to the odd channel interaction i 2 MD ^(l + 

x 2 MD ' ) P (T )pi 2 ■ 5(r 12 )p 12 , respectively Equation (jB13|) is available for the Skyrme interactions 
in which the LS currents are not ignored. 

We next show explicit form of the G™> factor in Eq. (jB10|) for typical interaction forms. 

1. 5 interaction 

Substituting f{2k\ 2 ) by 1, we obtain 

{h,k 2 ) =5 m . (B14) 
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2. Gauss interaction 

Because f(q) = (v^/aO 3 e~ {q ' 2 ^\ Eq. (lETTl) leads to 



+ m)! / /i 2 



For £ = and 1, we have 
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3. Yukawa interaction 

For the Yukawa interaction we use /(g) = 47r/ / u( / u 2 + g 2 ). Inserting it into Eq. 
we obtain for even £ 
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and for odd I 
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For £ = and 1, we have 
G< 0) (*i,k) 



vr /i 2 + (h + k 2 f 



fikikb /i 2 + (k 2 - /ci) s 



(B19) 
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3tt 



- 4fcifc 2 



(B21) 



Setting ki = k 2 = k F o and using the estimated level density at the Fermi momentum 
N = (2?t) 6 ■ 2/c F0 Mq/7t 2 , we define the usual Landau parameters 



f e = N Q F{ e) (k F0 , k F0 ) , f t = N FP (k F0 , k F0 ) , 
g t = NoF® (k F0 , k F0 ) , g> = N Q F$ (k F0 , k F0 ) . 



(B22) 



The second derivatives of £ at the saturation point are connected to the Landau param- 
eters. The following relations are verified, 



Ml 
M 



a x = 



6M r 



6M n * 



(B23) 
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